Disappearance of Transverse Flow in Central Collisions for Heavier Nuclei 
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Using the Quantum Molecular Dynamics model, we analyze the disappearance of flow in heavier 
colliding nuclei. A power law mass dependence (oc -^—) is obtained in all cases. Our results are in 

excellent agreement with experimental data which allows us to predict the balance energy for 238 U 
+ 238 U collision around 37-39 MeV/nucleon. 
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Thirty years ago it was predicted by Scheid and 
Greiner that in heavy ion reactions the nuclei will be 
compressed and heated and that this yields for non cen- 
tral reactions to in-plane flow (p* r ). More than a decade 
later, this conjecture was confirmed by the Plastic Ball 
group 0- In the following investigations it turned out 
that this in-plane flow carries information on the nuclear 
equation of state Q . If the nuclear equation of state is 
stiffcr, more compressional energy will be stored in semi- 
central reactions and, when released, more in-plane flow 
will be given to the nucleons. 

The maximal density which is reached in a reaction 
depends on the beam energy as well as on the system size. 
The lower the beam energy the less is the compression. 
At very low energies, the repulsive part of the nuclear 
equation of state, which appears at densities above the 
normal nuclear matter density, is not tested anymore and 
the nucleons feel only the attractive mean field. A typical 
example is the deep inelastic reactions in which the two 
nuclei rotate around a common center. This rotation 
creates in-plane flow as well but in opposite direction: 
Due to the common rotation the nucleons stick together 
for a while and will be emitted into the direction opposite 
to the impact parameter whereas the in-plane flow which 
is caused by compression will be in the direction of the 
impact parameter. 

There is a beam energy at which the in-plane flow dis- 
appears when changing from the direction into that op- 
posite to the impact parameter. It has been shown in 
the simulation of heavy ion reactions that this beam en- 
ergy called balance energy, E^ai, 0, IU El depends on 
the nucleon-nucleon (nn) cross-section in the medium as 
well as on the potential |(J. With the very recently 
measured Eb a i in 197 Au 4- 197 Au collisions Jj (earlier 
only estimated values were available @), a renewed in- 
terest has emerged in the field . In addition to the Au 
system, balance energies E ba i of 12 C + 12 C HH, 20 Ne + 
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, 40 Ar + 27 Al[l2, 40 Ar 
[3, 64 Zn + 27 A1 0, 40 Ar + 
58 Ni + 58 Ni 0,13, 64 Zn + 58 Ni 
93 Nb + 93 Nb||, 129 Xe + 118 Sn 
and 139 La + 139 La are also available. It is worth 
mentioning that most of the above studies were for the 
central collisions only. A few, however, also searched for 
the impact parameter dependence of the balance energy 



Ar + 27 A1 II 
13, 40 Ar + 51 V 

64 Zn + 48 Ti[nl 
Kr + 93 Nb HET 



Apart from the directed in-plane flow, differential as 
well as elliptic flow has also been predicted very recently 

G3 

The measurements of the balance energy over wide 
range of system sizes provide an excellent opportunity 
to pin down the role of the mass dependence, where 
only preliminary studies |?l.llO|| have been performed yet. 
These preliminary studies suggest a power law depen- 
dence ex A T of the balance energy on the mass number of 
the system. Interestingly, most of the theoretical stud- 
ies are done within the Boltzmann-Uehling-Uhlenbeck 
(BUU) model [!SSIl|lS|ll|ll|ll|l3|ll|lllM|23. 
Some attempts, however, also exist within the frame- 
work of Quantum Molecular Dynamics (QMD) model 
[l3l I22I 123L 01 ■ Heavy systems are rather rarely ana- 
lyzed in these approaches. 

Our present aim is therefore to study the mass depen- 
dence of the balance energy in heavy colliding nuclei and 
to predict for the first time the disappearance of the col- 
lective in-plane flow in central 238 U + 238 u collision. We 
shall show that the mass dependence of E ba i for heavier 
nuclei scales approximately more as rather than as 

A~i as has been suggested for light and medium collid- 
ing nuclei jlOj . The present study is made within the 
framework of QMD model which is described in detail in 

refs. n in in mm. 

In the QMD model, each nucleon propagates under 
the influence of mutual interactions. The propagation is 
governed by the classical equations of motion: 
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where H stands for the Hamiltonian which is given by: 
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Here vf kyrme , V? vk and Vf oul are, respectively, the 
Skyrme, Yukawa and Coulomb potentials. Momentum 
dependent interactions are not important at these low 
beam energies as Zhou et al. 21) have shown. It is worth 
mentioning that both the soft and hard equations of state 
have been employed in the literature. Following pl.l5l.IT2l 
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e he n mm mm Hi, we sha11 ais ° use a hard 

equation of state. For the nn cross-section we chose an 
isotropic and energy independent cross-section of 40 mb 
for the present analysis. This seems to us a reasonable 
choice in view of the fact that cross-sections based on 
G-Matrix calculations differ widely and because at this 
energy most of the collisions are Pauli blocked. Hence 
different cross-sections do not produce a large effect if 
they are not too small. 

Since we plan to study the heavier colliding nuclei, dif- 
ferent nn cross-sections should not have much effect |7J. 
Further, it has been shown in ref. 0] that the nucle- 
ons in the present energy domain (< 80 MeV/nucleon) 
collide with average er=55 mb. The isotropy of the 
cross-section also does not affect the reaction dynam- 
ics l28l . Similar assumptions were also made in refs. 
H El [H m m m H3- One should, however, keep 
in the mind that different nn cross-sections may affect 
the dynamics in lighter systems. 

Using the above description, we simulated the central 
collisions with system mass A (— At + Ap; At being the 
target mass, and Ap being the projectile mass) > f 75. In 
particular, we simulated 86 Kr + 93 Nb (b = 4.07 fm) E> 

93 Nb + 93 Nb ( b = 3 1Q4 fm ) [| j 129 Xe + 118 Sn ( b = q_ 3 

fm) |3, 139 La + 139 La (b = 3.549 fm) g, 197 Au + 197 Au 
(b = 2.5 fm) 7] and 238 U + 238 U (b = 0-3 fm) at incident 
energies between 30 MeV/nucleon and 80 MeV/nucleon 
at a step of 10 MeV. A straight line interpolation be- 
tween steps was used to calculate the energy of vanishing 
flow Eb a i. The reaction was followed till transverse flow 
saturates which is close to 300 fm/c for heavier collid- 
ing nuclei whereas it is < 200 fm/c for lighter colliding 
nuclei. 

In fig. 1, we display the average directed transverse 
momentum (p* r ) defined as: 

1 A 

<pf> = jEwjkw, (3) 

1=1 

where Y(i) and p x {i) are, respectively, the rapidity and 
the transverse momentum of the ith particle for the re- 
actions 93 Nb + 93 Nb, 139 La + 139 La, 197 Au + 197 Au and 
238 u _j_ 238 u yy e gee fo&t the collective in-plane flow for 
93 Nb + 93 Nb changes sign between 55-60 MeV/nucleon 
whereas it is already positive around 40 MeV/nucleon 
for 238 U + 238 u system. Further, the saturation time 
increases with the size of the colliding nuclei. The early 
onset of the flow in heavier colliding nuclei is due to the 
Coulomb forces that are much stronger in heavier systems 
compared to lighter nuclei. In addition, a large collision 
rate in heavier systems also contributes towards the early 
onset of the flow. 

In fig. 2, we display the energy of vanishing flow 
(Ebai) as a function of the combined mass of the sys- 
tem. As stated earlier, Eb a [ is extracted using a straight 
line interpolation between the calculated values of the in- 
plane flow. Here open squares represent our calculations 
whereas solid stars are the experimental findings. The 



dotted and dash-double-dotted lines are, respectively, the 
power law fits (cx A T ) to the theoretical values including 
238 U + 238 U in one case and excluding 238 U + 238 U in 
the other case. The fit to the experimental points is rep- 
resented by a solid line. All fits are obtained with \ 2 
minimization. The values of Eb a i (obtained with a stiff 
equation of state and 40 mb cross-section) are very close 
to the experimentally measured E^ai- The fit to the ex- 
perimental data yields r = —0.52451 ± 0.06261, whereas 
that to theory yields r = -0.53261 ±0.16373. Once 238 U 
+ 238 U is included, the r decreases to -0.50872±0. 10883. 
In other words, we observe a dependence in the 
Ebai- Similar dependence can also be obtained with a 
least square fit. Based on the above findings, we pre- 
dict Eb a i for the central 238 U + 238 U reaction around 
37-39 MeV/nucleon (according to power law fits, it is 
around 37-38 MeV/nucleon whereas QMD simulation 
predicts around 39 MeV/nucleon). It is worth mention- 
ing that most of the earlier mass dependence calcula- 
tions 0, E> E3 could not reproduce the experimentally 
extracted slopes 0, E ■ Zhou et al. [2lJ could repro- 
duce the slope, however their analysis was only done for 
lighter nuclei < 200. Our calculations can reproduce the 
experimentally extracted slope very closely, therefore, we 
can predict the energy of vanishing flow in the 238 U + 
U system. The large deviation of our r value from 
the standard value (« — -|) reflects the increasing im- 
portance of the Coulomb repulsion with the size of the 
system, as noted in ref. 01. There the r value is close 
to — | for masses < 200 EE3, whereas it increases to 
w —0.45 when heavier systems like 139 La + 139 La and 
197 Au + 197 Au are included. In the present analysis, we 
took only heavier nuclei (A > 175). Therefore the slope 
is steeper than the above cited values. If one also takes 
the lighter nuclei into consideration, our value also de- 
creases to —0.4. In other words, for lighter and medium 
nuclei, the balance energy E\, a \ emerges due to the in- 
terplay between the mean field and nucleon scattering. 
However, for heavier colliding nuclei, Coulomb interac- 
tion is as well an important factor. It is of interest to 
see the contributions of the mean field (that includes the 
Coulomb interaction) and nn collisions towards the trans- 
verse flow at the balance energy E^ai ■ Following ref. , 
we decomposed both these contributions in the simula- 
tions itself. At each time step during the reaction, the 
momentum transferred due to the two- body collisions 
and mutual mean field potential is calculated separately 
using equation (3) . The separation was done at each sim- 
ulated incident energy and a straight line interpolation 
was used. The decomposition is plotted in fig. 3 as a 
function of the total mass of the system. We find that 
the contribution of mean field towards transverse momen- 
tum is negative whereas it turns repulsive for the collision 
part. Both these contributions again obey a power law 
behavior with r = -0.70931 ± 0.24234. 

Summarizing, we present the disappearance of flow 
in heavier colliding nuclei with a prediction of balance 
energy for 238 U + 238 U around 37-39 MeV/nucleon. 
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Our calculations (with a stiff equation of state and 
o = 40 mb) are in a very close agreement with exper- 
imentally extracted values (jth — —0.53261 ± 0.16373; 
Texpt — —0.52451 ±0.06261). Both these findings suggest 
a power law mass dependence oc -^=. The contribution 
of the mean field towards flow at Eb a i is negative whereas 



it is positive for the collision part. Both contributions 
can be parameterized in terms of a power law. 
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Figure Captions 

FIG. 1. The time evolution of (p* r ) for four different 
reactions: (a) 93 Nb + 93 Nb, (b) 139 La + 139 La, (c) 197 Au 
+ 197 Au and (d) 238 U + 238 U. Here a stiff equation 
of state along with constant nn cross-section of 40 mb 
strength is used. 

FIG. 2. The E^ai as a function of the total mass of the 
system. Solid stars are the experimental data whereas 
open squares are the present theoretical results. The 



solid line is a \ 2 minimization fit of power law (oc A T ) 
for experimental data whereas dash-double-dotted line 
is a fit for the corresponding theoretical result. The 
theoretical fit that includes 238 U + 238 U reaction is 
represented by dotted line. 

FIG. 3. The decomposition of (p* r ) at E^ai into mu- 
tual mean field part and collision part as a function of 
the system size. The lines are the x 2 fit of power law 
oc c.A T . 
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Fig. 1 "Disappearance of Flow in Central Collisions" 
by Sood, Puri and Aichelin 
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Fig. 2 "Disappearance of Flow in Central Collisions" 
by Sood, Puri and Aichelin 
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Fig. 3 "Disappearance of Flow in Central Collisions" 
by Sood, Puri and Aichelin 



